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Abstract 

^ ■ By means of an updated renormalization method, we construct 

^*J , asymptotic expansions for unstable manifolds of hyperbolic fixed points 

(^ I in the double- well map and the dissipative Henon map, both of which 

^D ' exhibit the strong homoclinic chaos. In terms of the asymptotic ex- 

^1/ . pansion , a simple formulation is presented to give the first homoclinic 

Q>^ I point in the double- well map. Even a truncated expansion of the un- 

"^ ' stable manifold is shown to reproduce the well-known many-leaved 

J>-,' (fractal) structure of the strange attractor in the Henon map. 

V- 

o ■ 

=2 ■ 1 Introduction 

o ■ 

IJ • The bifurcation of separatrices and homoclinic or heteroclinic struc- 

tures are well known to lead to genesis of chaos in conservative dynamical 
systems |l[]. While numerical iterations of low-dimensional mappings easily 

5^ I provide these complicated structures, it is extremely difficult to derive them 

analytically. For nearly integrable systems, the difficulty may have been over- 
come by means of the asymptotic expansion beyond all orders supplemented 
with the Borel sum and the stokes phenomenon 00 0| p. However, this 
method can not be applied to the systems far from integrable such as a finite- 
time-discrete dynamical system with double-well potential (the double-well 
map) and the strange attractor in the dissipative Henon map, which will be 
analyzed in the present paper. 

Recently a novel method based on the perturbative renomalizaion group the- 
ory has been developed as an asymptotic singular perturbation technique P . 
The renormalization group (RG) method removes secular or divergent terms 



>$ 



from a perturbation series by renormalizing integral constants of lower order 
solutions. This RG method was reformulated on the basis of a naive renor- 
malization transformation and the Lie group [^ . The application of the RG 
method to some non-chaotic discrete systems was achieved in the framework 
of the envelope method 0. In this paper, we apply the reformulated RG 
method to chaotic discrete systems. Our approach is straightforward and 
much simpler than the original RG method [^ or the envelope method 
and we obtain asymptotic expansions of unstable and stable manifolds con- 
stituting homoclinic tangles in the double- well map and an unstable manifold 
included in the closure of the strange attractor of the dissipative Henon map. 

2 Double well map 

Let us analyze a symplectic map obtained by time-discretization of 
canonical equations for dynamical systems with a single freedom: 

Xn+l -Xn= Pn+1, Pn+1 - Pn = S^f{Xn)- (l) 

For concreteness, we choose a double well potential 
Then we get the second-order difference equation for x„ 

LXn = Xn+l - (2 + S^)Xn + Xn^l = -26^xJ, (2) 

where S'^ represents the time difference and is not a small parameter in 
general. This map is area-preserving and has a hyperbolic fixed point at 
{x,p) = (0,0). Although we concentrate on the system with a double- well 
potential in this section, the following analysis will hold for more general 
systems such as the standard map. For 6"^ ^ (i.e. in the continuous limit), 
the system becomes integrable and the phase space is occupied by regular 
orbits separated by a separatrix. If 5^ 7^ 0, the splitting of separatrix occurs 
and unstable and stable manifolds of the hyperbolic fixed point cross each 
other at an infinite number of points (homoclinic points) that accumulate 
to the hyperbolic fixed point. This complex homoclinic structure (homo- 
clinic tangles) is the typical chaos generated by the Birkhoff-Smale's horse 
shoe mechanism. When 6"^ is small enough, the splitting of separatrix is 
very small and an asymptotic expansion of unstable and stable manifolds 



was obtained by means of the singular perturbation method of beyond all 
order 0. Here, we are interesting in the case that 5"^ ~ C^(l) or larger. 
Let us construct a formal series solution x^ near the hyperbolic fixed point 
(a;,p) = (0,0) along the unstable manifold, that is, lim„^_ooa;5^ = 0: (For 
brevity, x^ will be noted as Xn in the following.) 

^ _ ^(1) I ^(2) I (3) I (4) J (o\ 

Here 5^ is not necessarily small and \xn \ > \xn \ > ■■■. Then, naive per- 
turbed equations for Xn are written as 

Lx« = 0, Lx^^^=0, 

r^(3) _ _OX2^(l)3 r^(4) _ _^r2 (1)2 (2) 



The leading-order solution satisfying the condition lim„^_ooa;„ = is given 

by 

xW = AK\ (4) 

where A is an integral constant, 6 = 2sinh(/c/2) > and K = expk > 1. 
The series solution (||) is obtained as 



Xr. 



AK"" (l - C2A^K^'' + CiA^K^"" - cqA^K^'' + ■■■), (5) 



where C2, C4 and cg are constants defined as 



36^ 1 _3_ 



ce = ——-(—- + —-), Dj- = coshj/c-coshA; ij ^ N). 



Now we introduce a renormalized variable A{n) so that all the higher order 
terms in the series solution (|^) are renormalized into the integral constant 
A, that is, 

Xn = i(n)ir^ (6) 

where 

A{n) = a{i- C2A^K^'' + c^A^K^'' - cqA^K^'' + ■■■), (7) 

which is called a renormalizaion transformation A ^ A |p. In order to 
derive a difference equation (a renormalizaion equation) for the renormalized 



variable A{n), we express a difference A{n + 1) — A{n) in terms of A. From 
d^), we have 

A{n+1)-A{n) = -C2{K^ - l)K'^''A{nf 

+Ci{K^ - l)K^''A{n)'' + ■■■. (8) 

By solving (0) with respect to A iteratively, A can be expressed in terms of 
A and n as 

A = A{n){l + c^AinfK^'' + (Sc^ - C4)i(n)^ir^" + ■■■}. 

Thus we can replace A in (|]) by A and obtain a renormalizaion equation: 

i(n + l)-i(n) = -C2{K^-l)K'^''A{nf 

+{ci{K^ - 1) - 3c2(K2 - l)}K^"i(n)^ 

-{C6(ir6 - 1) - 5c2C4(ir^ - 1) 

+3c2iK^ - l)(4c2 - C4)}ir6"i(n)^ + ■ ■ ■ , (9) 

which is a non-autonomous difference equation. In terms of the original 
variable x„ (see (^), the renormalizaion equation (|^) is transformed into an 
autonomous one. 

Xn+l = KXn- KC2{K^ -l)xl 

+K{c4{K^ - 1) - 3cl{K^ - l)}x^ 

-ir{c6(/C6 - 1) - 5C2C4(ir4 - 1) 

+ 3c2{K^-l){4:cl-C4)}xl + --- 

= giM, (10) 

which gives the first branch gi of a return map on the unstable manifold: 
Xn+i = gj{xn), where gj is a multi-valued function of x„ and j(= 1,2, ■ ■ •) 
denotes a branch number. When K ^ 1 (strong chaos) , it is easy to see 
that 

g,{x) = K[f{x) + 0{K-% (11) 

in general. If {K — 1) ^ 1 (weak chaos), we have, 

g,{x) =x + {K- l)x(l - x^y/^ + 0{{K - 1)^). (12) 

for the double well map for i^ = 2.1, an asymptotic expansion oi gi truncated 
at xj^ is depicted in Fig .1. Although the truncated expression agrees well 



with the exact numerical result for x„ < 0.75, it deviates considerably from 
the exact result near a singular point or a branch point (x„ ~ 1.1) . In order 
to avoid this considerable discrepancy near a branch point, we restrict the 
domain of gi so that the return map gi is reversible, that is, Xn ^ 0.75. The 
reversible branch thus obtained is denoted by gi . Then, from @, we have 
the following functional map: gj -^ cjj+i 

gj+l{Xn) = 2Xn - gJ^{Xn) + 5'^f{Xn), (13) 

where the domain of Xn is chosen for each j so that a new branch ^j+i is 
reversible. Using the functional map (|13|), we can construct each reversible 
branch of the return map (jj from the first branch gi or gi step by step. As 
shown in the figures 2 and 3, the result agrees well with the exact one even 
if a truncated expansion of gi up to x^ is used for gi. Once the return map 
cjj on the unstable manifold is constructed, the unstable manifold (x",p") is 
written as 

(14) 



p" 


= x^-9] 


-\x 


ecause pl = xl- xl_^ = x^- 


-97\x-). 




Using the symmetry of the 


map (|l|): 




x' = 


x-p, 


V- 



-p, (15) 

the stable manifold (x*,p*) is constructed as follows. The symmetry (|T5p 
gives 

x" = x'-p^ p'' = -p'. (16) 

From ([T6|) and (0), we obtain the stable manifold in terms of gj 

f = x'-~g,{x'). (17) 

Therefore, the x-coordinate of a homoclinic point {x^,p^) is given by 

x'' = ~gfg,{x''). (18) 

Since the first intersection point {x*,p*) is located at p* = Q on the first 
branch gi or ^2 in the present case, x* is a fixed point of the first branch of 
the return map. 

x* = g,{x*) = ~g2{x*). (19) 

For /r > 1 and (iT - 1) < 1, (0), (0) and (0) yield x* ^ ±l/v^ and 
X* ~ ±1 respectively. 



3 Henon map 

Here, we construct an unstable manifold in the strange attractor of a 
dissipative dynamical system known as the Henon map: 

Xn+i = 1-axl + yn, (20) 

y„+i = bxn, (21) 

where a = 1.4 and b = 0.3. This map has the hyperbolic fixed points: 

, _ b-l±{{b-lf + Aay/^ 

2a 
y* = bx*. 

Let us consider the unstable manifold of one of the fixed points: x* = [6 — 1 + 
{(6—1) + 4a}^/^]/2a, y* = bx*. This fixed point is known to be included in 
the closure of the strange attractor of the Henon map. Setting Xn = Xn — x* 
and eliminating the y-variable, we obtain the second-order difference equation 
from (EOD and (M)- 



x„+i + 2ax*Xn - bxn-i = -ax\. (22) 

Near the origin of £, we obtain the following formal series solution which 
vanishes as n ^ — oo. 

x„ = K''[A- 02^2 TT" + csA^iT^" - c4A^ir=^" + ■■■], 

where A is an arbitrary integral constant; K = —ax* + {(ax*)^ + b}^^"^ ~ 
—1.923739 is one of the eigenvalues of the linearized map of (^) and 

aK^ 2a^K^ 

C2 = TTTFT^^ C3 



a^K^ 1 AK 

D{K^) = K^^ + 2ax*K^ - b. 



Following the same procedure as in the previous section, we renormalize A 

as 

A{n) = A-C2A^K'' + C3A^K^''-C4A^K^'' + ---, 

Xn = i(n)K", (23) 

6 



and obtain a similar non-autonomous renormalization equation of A as (^. 
Replacing A{n) by x„ in the renormalization equation of A, we have 

Xn+l = KXn - C2xlK{K - 1) 

+xlK{K-l)[c,{K+l)-2cl] 

+xiK{K - 1)[-C4(ir2 + K + l) + 3c3C2(ir + 1) 

-C2(5c2 - cl)] 
+xlK{K - l)[c^{K^ + K^ + K+l)- Ac2Ci{K^ + K + 1) 

+3C3(3C^ - C3)(ir + 1) - 2c2{lcl - 6C2C3 + C4)] + • • • 

= giixn), (24) 

which gives the first branch gi of a return map on the unstable manifold: 
Xn+i = gj{xn), where gj is a multi-valued function of Xn and j(= 1,2, ■ ■ ■) 
denotes a branch number. An asymptotic expansion of gi truncated at x^ 
is depicted in Fig. 4, where the truncated expression deviates considerably 
from the exact numerical result near x„ ~ 0.7 and £„ < — 1 . In order 
to avoid this considerable discrepancy, we restrict the domain of gi{x) to 
—0.7 > X > 0, where the return map gi is reversible. The reversible branch 
in this way obtained is denoted by gi . Then, from (p2D, we have the following 
functional map for neighboring reversible branches: cjj -^ gj+i 

gj+i{xn) = -ax^ - 2ax*Xn + bgj^{xn), (25) 

where the domain of x„ is appropriately chosen for each j. Using the func- 
tional map (P5|), we construct each reversible branch of the return map gj 
from the first branch gi (0) step by step. Even when the asymptotic ex- 
pansion of gi is truncated at x^, the result agrees well with the exact one 
as shown in Fig. 5 and Fig. 6. In Fig. 7, we see how well each branch of 
the return map constructed from an initial truncated map gi recovers the 
well-known many-leaved (fractal) structure of the strange attractor near the 
fixed point . 

4 Conclusion 

We construct asymptotic expansions for unstable manifolds of hyperbolic 
fixed points of the double-well map and the dissipative Henon map, both of 
which exhibit the homoclinic chaos. Since the dimension of the manifold is 



one, the dynamics on the manifolds is described by a return map with an 
infinite number of branches. An asymptotic expansion of the first branch of 
the return map is obtained by means of an updated renormahzation method 
which is an extension of the RG method. We exphcitly give a functional 
map, through which the other branches are obtained consecutively from the 
first branch. Thus a global asymptotic form of the unstable manifold is 
constructed to give the homoclinic chaos in the double-well map and the 
strange attractor in the dissipative Henon map. Using a truncated expansion 
of the first branch of the return map, we calculate an approximate unstable 
manifold, which agrees well with the exact numerical result. A fixed point of 
the first branch of the return map is found to be the first homoclinic point 
in the double well map, and the approximate unstable manifold reproduces 
the well-known many-leaved (fractal) structure of the strange attractor in 
the Henon map. 
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Figure 1: The first branch of the return map on an unstable manifold for 
the double well map. While dots ( • ) denote an exact solution calculated by 
means of numerical iterations, pluses (+) denote gi obtained by the renor- 
malization method. 
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Figure 2: The exact return map on an unstable manifold for the double well 
map. 
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Figure 3: The return map plotted by means of matching the set of Qj. 
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Figure 4: The first branch of the return map on an unstable manifold for 
the double well map. While dots ( ■ ) denote an exact solution calculated by 
means of numerical iterations, pluses (+) denote g\ obtained by the renor- 
malization method. 
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Figure 5: The exact return map of the Henon map's unstable manifold is 
plotted by numerical iterations. 
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Figure 6: The return map plotted by means of the matching the some of 
the (jj . 
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Figure 7: The return map of the Henon map's unstable manifold near the 
hyperbolic fixed point (a:„, i:„+i) = (0, 0) where the fractal structure appears. 
Pluses (+) express the exact map and line express the some of the Qj. 
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